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ABSTRACT 

Aims. We present the analysis of the [a/Fe] abundance ratios for a large number of stars at several locations in the Milky Way bulge 
with the aim of constraining its formation scenario. 

Methods. We obtained FLAMES-GIRAFFE spectra (R=22,500) at the ESO Very Large Telescope for 650 bulge red giant branch 
(RGB) stars and performed spectral synthesis to measure Mg, Ca, Ti, and Si abundances. This sample is composed of 474 giant stars 
observed in 3 fields along the minor axis of the Galactic bulge and at latitudes b=-4°, b=-6°, b=-12°. Another 176 stars belong to a 
field containing the globular cluster NGC 6553, located at b=-3° and 5° away from the other three fields along the major axis. Stellar 
parameters and metallicities for these stars were presented in Zoccali et al. (2008). We have also re-derived stellar parameters and 
abundances for the sample of thick and thin disk red giants analyzed in Alves-Brito et al. (2010). Therefore using a homogeneous 
abundance database for the bulge, thick and thin disk, we have performed a differential analysis minimizing systematic errors, to 
compare the formation scenarios of these Galactic components. 

Results. Our results confirm, with large number statistics, the chemical similarity between the Galactic bulge and thick disk, which 
are both enhanced in alpha elements when compared to the thin disk. In the same context, we analyze [ayFe] vs. [Fe/H] trends across 
different bulge regions. The most metal rich stars, showing low [ayFe] ratios at b=-4° disappear at higher Galactic latitudes in agree- 
ment with the observed metallicity gradient in the bulge. Metal-poor stars ([Fe/H]< -0.2) show a remarkable homogeneity at different 
bulge locations. 

Conclusions. We have obtained further constrains for the formation scenario of the Galactic bulge. A metal-poor component chemi- 
cally indistinguishable from the thick disk hints for a fast and early formation for both the bulge and the thick disk. Such a component 
shows no variation, neither in abundances nor kinematics, among different bulge regions. A metal-rich component showing low [a/Fe] 
similar to those of the thin disk disappears at larger latitudes. This allows us to trace a component formed through fast early mergers 
(classical bulge) and a disk/bar component formed on a more extended timescale. 

Key words. Stars: abundances, late-type - Galaxy: abundances, bulge, formation 



1. Introduction 

Nearly 25% of the visible light emitted from stars in the lo- 
cal universe c omes from bulges embedded in galactic disks 
dFukugita et al. Therefore, a fundamental step in the un- 

derstanding of galaxy formation and evolution is to obtain an ac- 
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curate characterization of these structures. The picture seems to 
be complicated as recent observations suggest that bulges may 
be an in-homogeneous class of ob jects, each of them holding 
a different formation history (e.g. lAthanas soula 2005). In this 
context, The Milky Way bulge is a unique laboratory in which 
we can study in detail the stellar populations as well as the 
chemical composition an d kinematics on a star-by-star basis (see 
iMinniti & Zoccalill2008l for a review). 

In our galaxy, the consensus is that the bulk of the bulge is 
old (t~12 Gyr) with a metallicity distribution that extends from 
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[Fe/H"|=-1.5 to rFe/H1=0.5 , with a peak at solar m e tallicity (e.g. 
i 



Zocc ali et all 120031 120081: IClarkson et aD 120081: iBrown et al l 
2010, and references therein) and that the formation of the 
bulge was fast as evidenced by the enhanced alpha elements 
in bulge stars dMcWilliam & RichlJT994t iRich & Origlial 120051; 



200 71: ILecureur et alj 
Rvde etalJ I2009L boiOt). 



Zoccali et al. 2006; Fulbright et al 

20071: IMelendez et all 120081: iRvde et alJ 12009. 2010). A verti- 
cal metallicity gra dient is also observed (Mi nniti et al.l Il995t 
IZoccali et al]|2008l) al though it appears to be absent in the inner 
bulge regions (b<4°) dRich et alJ f2007). Therefore, the conclu- 
sions so far regarding age, metallicity and chemical enrichment 
in the bulge point to a formation scenario similar to that of early- 
type galaxies through one or more starburst events triggered by 
mergers in the early ages of the universe dWeinzirl et al.l f2009; 
Hopkjnsetal|2010) - in other words, as a classical bulge. 

On the other hand, Kormendv ( 1993) presented observations 
of a large number of bulges with properties resembling more 
disk-like structur es rather than the spheroidal c oncept of clas- 
sical bulges (see iKormendv & Kennicutl [2004. for a review). 
These structures, named pseudo-bulges or disky bulges, tend to 
show younger stellar populations, to be rotation supported, and 
to have surface brightness profiles similar to those of disks and 
less concentrated than those of classical bulges. They are ex- 
pected to be formed from gas inflows driven by disk instabil- 
ities such as the presence of a central bar, which are common 
phenomena in spiral galaxies. Galaxies, when observed edge- 
on, often show a central component which s wells out of the disk 
with a boxy, peanu t, or even an X shape (e.g. lBureau et al.l2 006; 
Patsi s et al.ll2002l) . Several numerical simulations have shown 
that the inner parts of a bar might evolve and buckle off the plane 
of the disk showing these particular sha pes (Athanassoula 2005; 
iDebattista & Williamsll2004l;lDebattista et alj|2006l) . 

The above scenarios for bulge formation have historically 
emerged from the observation of bulges in nearby galaxies. 
However, a new paradigm has recently emerged, mostly driven 
by the observations of galaxies at redshift ~ 2, i.e., at a look- 
back time comparable to the age of bulge and thick disk stars. 
Such galaxies may offer examples of how our own Milky Way 
was looking like in its early days, prior and during its bulge for- 
mation. At z ~ 2 a widespread population exists of large, rotat- 
ing disk galaxies, with much higher g as fractions and velocity 
dispersion compared to local spira ls (Genz eTet al.ll2006l [2 008: 
Forster Schreib er et al.l 120091: iTacconi et al.l 1201 Ol iDaddi et al.l 



20101) . Such gas-rich galaxies are prone to disk instabilities that 
can result in bulge formation over timescales of a few 10 s yr, 
much shorter than those t ypically ascribed to secular instabilities 
in local disk galaxies dlmmeli et al] 12004 ICarollo et all 120071: 
lElmeg reen et al. 200H). Such rapid formation, partly due to the 
quasi-exponential mass growth experienced by these galaxies, 
quite naturally leads to a n q-element enhancement in th e re- 
sulting stellar populations dRenzinill20"09"l: iPeng et al.ll2010l) . that 
would apply both to the disk and bulge stars formed at these 
early cosmic epochs. 

While the existence of a metalli city gradient in the bulge, 
as shown by IZoccali et alJ (120081) strongly argues agains t 
the pseudo-bulge hypothesis of IKormendv & KennicutJ ( 12004) . 
some recent chemical abundance studies and morphological 
signatures have provided supporting evidence for such a sce- 
nario. Using a sample of hig h resolution spectra for disk and 
bulge giants, Mel endez et al J (l2008h find a similarity between 
the [O/Fe] ratio of the bulge and the thick di sk, both enhanced 
when compared to the thin disk. Later on, lAlves-Brito et al.l 
(l20Toh reach the same conc lusions when a n alyzin g other a- 
elements from the sample of Mel endez et al] d2008l) . These re- 



sults are in disagreement with previous results which show a 
bulge more enhanced in alpha elements than the thick disk. Since 
those previous studies were don e comparing bulge giants to local 
disk dwarfs from other work s (Bensbv et al. 2005; Reddv ~et~aT] 
2006), IMelendez etail (2008) claim that differences between 
thick disk and bulge stars might be explained by systematic er- 
rors when comparing dwarfs to giants or relative calibrations be- 
tween those works. If the similarity between the bulge and thick 
disk is confirmed with a larger and homogeneous sample, this 
would allow for a common origin for these two galactic compo- 
nents, during the early life of the Galaxy. However, it is impor- 
tant to recall that the origin of the thick disk i s also unknown, 
thoug h several hypothesis have been explored (Di Matteo et al. 
1201 ll and references therein), and most of the results are based 
on stars from the solar neighbourhoo d. Progress on this point 
has been achieved for the first time bv lBensbv et al.l (12010) who 
provided a sample of inner disk stars which show no significant 
chemical differences with respect to local disk giants. 

In terms of morphology, the evidences for a barred struc- 
ture in the Milky Way are numerous. Red clump stars distri- 
butions from photometric surveys have established the exis- 
tence of a bar inclined with respect to the Sun-Galactic cen- 
ter line of sight with its near end towards positive longitudes. 
Observations of the bulge are usually modeled with a triaxial 
bulge with a positi on angle of ~ 15-30 degrees and nearly 2.5 
kpc in length (e .g. Babusiaux & Gilmore 2005; R attenburv et al] 
l2007l) . lHoward et al.l(l2009l) analyzed the kinematics of the bulge 
concluding that the dynamical signature of the stars at high 
latitudes (b=-8°) is the same as that in fields closer to the 
galactic plane (b=-4°), therefore showing evidence for cylindri- 
cal rotation as expected for boxy bulges. Their data lack evi- 
dence indicative of an accretion origin of bulge in contrast with 
the claims of tidal streams in kinematics and metallicity from 
Rang wala & W illiams (2009) . Moreover thei r cylindrical rota- 
tion signal was modelled bv lShen et al] (1201 Ol) without any need 
for a classical component. 

There has also been evidence, from both observations and 
models, that different types of bulges could coexist in the same 
galaxy dSamland & Gerhard! [20031: iNakasato & Nomotd l200l 
IKormendv & Kennicutlll20041 : lPeletier et al]l2007l) . If this is the 
case, we might have a very complex structure in the bulge, with 
two or more components partially overlapping. 

Evidence for a dual nature of our bulge has most recently 
been invoked based on chemical abundances and kinematic stud- 
ies. Indication for bimodality in the metallicity distribution has 
been found from red clump stars in Baade's Window (Hill et al. 
201 1 submitted ). The two components seem to have also differ- 
ent kinematics (Babusia ux et alJ l2010) which is consistent with 
different structures present in the bulge: a metal-rich with bar- 
like kinematics and a metal-poor one showing kinematics of a 
spheroidal component. 

The picture is further complicated by the recent finding of a 
double peaked red clump along the minor axis, at both positive 
and negative latitudes, but only for | b |> 6°. The double clump 
is clearly visible in seve ral color mag nitude diagrams avail- 
able in the literature (e.g. [Zoccali et al.l 12003 ] : [Rattenburv et alj 
I2007L lMcWilliam&Zoccalill2010HNataf etalJ l2010h . Such a 
feature has been interpreted as two components at different dis- 
tances and is likely to be th e outcome of an X-shaped bulge 
dMcWilliam & Zoccalill2010l) . 

All these observational evidences have suggested a very 
complex structure of the bulge, with two or more distinct compo- 
nents partially overlapping. Signatures of the different formation 
histories for these components, should still remain in bulge stars. 
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Table 1. Galactic coordinates of bulge fields. Extinction values 
and number of stars analyzed in each field are also listed. 



N 


Name 


r 


b° 


E(B-V) 


stars 


1 


Baade's Window 


1.14 


-4.18 


0.55 


194 


2 


b = -6° 


0.21 


-6.02 


0.48 


194 


3 


Blanco field 


0.00 


-12.0 


0.20 


86 


4 


NGC 6553 


5.25 


-3.02 


0.70 


176 



In this paper we analyze the a-element distribution of gi- 
ant stars in three fields along the minor axis and another at a 
larger galactic longitude. We use them to homogeneously com- 
pare with disk giant stars abundances and to investigate the pres- 
ence of gradients within the bulge. 



2. The sample 

The bulge s pectra discuss e d here belong to the same sample 
presented in IZoccali et alJ d2008, hereafter Z08) as part of a 
project aimed to provide a spectroscopic characterization of the 
Galactic bulge. The sample consists of spectra of 650 K gi- 
ant stars obtained w ith the FLAMES -GIRAFFE spectrograph 
dPasquini et al.ll2003[) using the HR13, HR14 and HR15 setups 
(6100-7000A) at a resolution of R~ 20,000. For the present anal- 
ysis we have only used the HR13 setup (R=22,500) in which we 
have lines for all the analyzed elements. As described in Z08, 
the S/N of our sample ranges from 40 to 90. Details on the fields 
analyzed here are presented in Table [TJ In the last column of 
Table QJ we list the number of stars for which we were able to 
measure abundances for Ca, Mg, Ti and Si. These are only those 
spectra that had sufficient S/N (>50) for all the elements to be 
measured and did not present very large FWHM which could 
affect the results due to line blending. The sample also contains 
target stars which belong to bulge globular clusters: NGC 6522 
in Baade's Window (7 stars), NGC 6558 in the field at b=-6° (3 
stars) and NGC 6553 in the eponymous field (29 stars). These 
stars were clasified as cluster members by Z08 based on the fol- 
lowing criteria: [Fe/H] within 0.2 dex and radial velocity within 
+ 10 km/s of the mean for the cluster, and a location within 2 
arcmin from the cluster center. 

Magnitudes and astrometry for this sample come from the 
OGLE catalog for Baade's Window and from WFI images for 
the field at b=-6°, b=-12° and NGC6553 (see Z08 for a detailed 
description). Selected targets are nearly 1 mag above the red 
clump as shown in FigQJ 

The sample for thick and thin disk RGB stars in the lo- 
cal neighborho o d con sists of the same spe c tra an alyzed in 
iMelendez etaf] d2008l) and lAlves-Brito et all d2010l hereafter 
AB10). Stars belong to a range of effective temperatures of 3800 
< T e ff < 5000 K and metallicities -1.5 < [Fe/H] < 0.5, mak- 
ing this sample ideal to compare with our bulge sample. Spectra 
were taken with different spectrographs, MIKE (R~65,000) 
mounted on the Clay 6.5m Magellan telescope, the 2dcoude 
(R~60,000) spectrograph on the 2.7m Harlan J. Smith tele- 
scope at McDonald Observatory, HIRES (R~ 100,000) on the 
Keck I 10m telescope and complemented with spectra from the 
ELODIE archive (R~42,000). Population membership and data 
reduction is described in AB10. Population membership for the 
bulge sample is described in Z08. 
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Fig. 1. Color magnitude diagrams for the four fields: Baade's 
window (upper left), b=-6° [upper right), b=-12° (lower left) 
and NGC6553 (lower right). Spectroscopic targets in each field 
are marked as large filled circles. Magnitude s were obtained 
from the OGLE catalogue (lUdalski et al.ll2002l) and the Z08 cat- 
alog obtained from WFI images 



3. Stellar parameters 

For bulge stars we have adopted the same stellar parameters, ef- 
fective temperature (T e ff), surface gravity (logg), microturbu- 
lence velocity (£) and metallicity [Fe/H] presented in Z08 and we 
refer the reader to that work for a detailed description of the pro- 
cedure. In order to perform a differential analysis when compar- 
ing bulge and disk stars, we followed the same procedure from 
Z08 to determine parameters of the disk sample. We have used 
the same iron linelist from Z08 based on calibrated log gf values 
to match the observed spectra of the metal-rich giant //-Leonis. 
Only small difference is that as first input value for T e ff we 
us ed photometric value s obtained from th e calibration provided 
bvlAlonso et alJ d 19991) as opposed to the iRamfrez "& Melendez 
d2005h calibration adopted by Z08. T e jf values were then refined 
spectroscopically by imposing excitation equilibrium between 
Fel lines therefore removing any dependance on the adopted 
calibration. Microturbulence velocity £ was calculated by re- 
quiring the same [Fe/H] for all equivalent widths of Fel lines. 
Surface gravity was determined photometrically using the fol- 
lowing equation: 



log (g) = log 



M 

So ) + log ( — 1 + 0.4 (M M ,» 



M So/ , ) + 41og 



(1) 
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Table 2. Atomic linelist for Ca, Mg, Ti and Si used in this 
work. Also listed are the excitation potential ix ex ) and oscilla- 
tor strength (log gf) for each analysed line. 



/1(A) 


Element 


Xex 


logg/ 


6318.71 


Mgl 


5.110 


-2.000 


6319.23 


Mgl 


5.110 


-2.240 


6319.49 


Mgl 


5.110 


-2.680 


6166.43 


Cal 


2.521 


-1.142 


6169.04 


Cal 


2.523 


-0.797 


6169.56 


Cal 


2.526 


-0.478 


6312.23 


Til 


1.460 


-1.552 


6142.49 


Sil 


5.620 


-1.500 



Table 3. Abundances of the Sun, Arctums and //-Leonis, ob- 
tained by comparison to synthetic sp e ctra, a s well as reference 
solar abundances from lAsplund et ali ([2009). The last two rows 
show the zero point abundances adopted by AB10 and our mean 
abundances measured for those same disk giants to which they 
refer. 



Stars 


A(Ca) 


A(Mg) 


A(Ti) 


A(Si) 


Sun 


6.39 


7.60 


4.95 


7.51 


Sun (Asplund et al) 


6.34 


7.60 


4.95 


7.51 


Arcturus 


5.93 


7.30 


4.55 


7.17 


/j-Leonis 


6.65 


8.02 


5.31 


7.85 


AB10ZP 


6.32 


7.65 


4.83 


7.60 


This work 


6.49 


7.65 


5.09 


7.62 



where for the Sun we have used the same values used by Zoccali 
et al. (2008): bolometric magnitude Mb i,g = 4.72 mag, effective 
temperature T e jf = 5770 K and surface gravity log (g Q ) 
4.44 dex. We have adopted the masses from AB10 which are 
based on Padova isochrones. Absolute visual magnitudes were 
obtained assuming distances from Hipparcos parallaxes. Finally, 
bolometric mag nitudes were calcula ted by applying bolometric 
corrections from lAlonso et al. I (fl999l 



4. Abundance analysis 

Mg, Ca, Ti and Si abundances were determined by com- 
paring th e observed sp ectra with synthetic ones created with 
MOOG (ISnedenlll973l) . MOOG is a FORTRAN code that per- 
forms a spectral line analysis and spectrum synthesis assum- 
ing local thermal equ ilibrium. MARCS model atmospheres 
dGustafsson et al. 2008) were used for our analysis. These mod- 
els were created using the stellar parameters provided in Z08 for 
bulge stars, while for disk stars we used the stellar parameters 
determined as described in Section 3. We have used the same 
atomic line list from the analysis of bulge stars in lLecureur et al.l 
d2007l hereafter L07). Unfortunately, the lower resolution of 
GIRAFFE compared to that of the UVES spectra of L07 pre- 
vents us to obtain accurate oxygen abundance measurements 
from the available OI line at 6300 A. 

The broad (~5 A) autoionization Ca I line at 6318.1 A affects 
the continuum of the Mg triplet as discussed in L07. Therefore, 
Ca abundances were determined first in our procedure and then 
this value was given as an input to determine Mg abundances. 
The CN line list from the same work was adopted and calibrated 
carefully as CN lines are known t o be impo rtant for cool giant 
stars. The TiO molecular line list (|Pleil998) was also included. 
Reference solar abundances are from Asplun d et al.l (12009). 

Although this is a differential analysis, in the sense that 
we determine the stellar parameters and perform the spectral 
analysis in the exact same way for all stars, it is important to 
provide abundances of some reference targets using our pro- 
cedure, therefore establishing our "zero point" in order to be 
able to contrast our results with the values obtained in other 
works. For this reason we have also adopted the calibration 
procedure based on reproducing the observed spectrum of the 
Sun, Arcturus and jif-Leonis. We have used the optical spectra 
of these stars that were analyzed in L07 for the same purpose. 
//-Leonis spectrum was taken at the Canada-France-Hawaii with 
the ESPanDOnS spectropolarimeter at a resolution of 80,000. 
The Arcturus spectrum with a resolution of 120,000 comes from 




6311.6 6312.0 6312.4 6312.8 
Wavelength (A) 



6141.6 6142.0 6142.4 6142. £ 
Wovelength (A) 



Fig. 2. Calibration of the synthesis regions for the Sun, Arcturus 
and yU-Leonis. Observed spectra convolved to the GIRAFFE res- 
olution are shown as black circles and the red solid line shows 
the synthetic spectra using the stellar parameters from L07 and 
abundances in Table [3] 



the UVES database (Ba gnulo et_ alj2003h . Finally for the Sun we 
used a UVES optical spectrurrQ. These spectra were convolved 
to the GIRAFFE resolution and compared with synthetic spectra 
in a region around the Mg, Ca, Ti, and Si lines analyzed in this 
work. Table [2] shows the wavelength, excitation potential, and 
log gf of these lines. Model atmospheres were created adopting 
the same stellar parameters used in L07, [Fe/H]=0.3, T e fr =4540, 
£=1.3, logg=2.3 for^-Leonis, [Fe/H]=-0.52, T eff =4300, £=1.5, 
log g= 1.5 for Arcturus and T eff =5770, £=0.9, logg=4.42 for the 
Sun. Figure [2] shows the comparison between the synthetic and 
observed spectra for these stars with the abundances listed in 
Table [3] For the Sun, we only required a 0.05 dex modifica- 
tion in Ca abundanc e with respect to our reference values of 
lAsplund et ail (120091) also listed in Table |3] 

Abundances for our target stars were determined by an iter- 
ative process in which the synthetic spectrum is compared to the 
observed one modifying the abundance in each step of the pro- 
cess until reaching the best fitting value. Given the large number 
of stars in our sample, we have developed a code which calcu- 
lates the best fitting abundance by x 1 minimization in a semi- 
automatic way. In the first step, radial velocities and a first guess 



(http://www.eso.Org/observing/dfo/quality/UVES/-pipeline/solar_spectrum.l 
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6311.6 6312.0 6312.4 6312.8 
Wavelength (A) 



6'41.6 6142.0 6142.4 6142.8 
Wavelength (A) 



Fig. 3. Example of the final output of the synthesis procedure 
in order to derive the abundances of Mg (upper left), Ca (upper 
right), Ti (lower left) and Si (lower right). Observed spectrum for 
target star 212 175c6 is plotted as black dots and the red solid line 
shows the synthetic spectrum using the best fitting abundance as 
obtained from our procedure. 



for the FWH M of each star is d etermined using the DAOSPEC 
code dStetson & Pancinoll2008l) . The placement of the contin- 
uum is carried out manually for each star, by direct comparison 
between observed and synthetic spectrum in a region of 10 A 
around the line of interest. Additionally, inspection of nearby 
lines are used to modify the FWHM from the value obtained 
from DAOSPEC, when the latter clearly fails to reproduce the 
spectral features. Later on, the code automatically calculates the 
best fitting abundance by minimization of a x 2 value which is 
obtained in the following way: 



W(/l) 



(2) 



where F sy „(A) is the flux of the synthesis at each wavelength, 
Fobs(X) is the observed flux previously normalized and W(A) is 
a weight value calculated as the difference between a normal 
synthesis and another without the line of the element to be fitted. 
This allows us to give more importance to the core of the line un- 
der analysis. To create this weight we h ave used n Leonis param - 
eters and abundances as a reference dGratton & Sneden l ll990l) . 
Figure [3] shows an example of the final output from the code for 
each analyzed line. 

5. Error analysis 

5. 1 . Error due to stellar parameters uncertainties 

We have carefully analyzed the effect of the different sources 
of error in our procedure. When obtaining stellar abundances 
the errors are commonly dominated by the uncertainties in the 
stellar parameters. Their influence can be determined by vary- 
ing the adopted parameters by their respective uncertainty and 
then re-calculating the abundance value which best fits the ob- 
served spectra. We have carried out this analysis in a subsam- 
ple of stars which covers the whole range of metallicities and 
effective temperatures of our sample. We have used the uncer- 
tainties in stellar parameters presented in Z08 of 200K for T e ff , 
0.3 dex in logg, 0.2 dex in £ and 0.1 dex in the final [Fe/H] 
value. Figure [4] shows the variation in determined abundances 
for each element when varying the stellar parameters by their 



respective uncertainty. The final errors (r r ,cr [Fe) / H] ,cr logf ,,crf asso- 
ciated to each parameter were calculated as the average between 
the mean changes in abundance when uncertainties were added 
(<t (+) ) and subtracted (cr ( _)). 

Figure [4] shows that errors in the abundance determination 
are strongly dominated by the uncertainties in the effective tem- 
perature determination. Given the adopted uncertainty of 200 
K in temperature, Mg seems to be the least sensitive element 
with an associated error of <x Ar =0.05 dex. Ca and Si show an 
associated error of 0.07 and 0.10 dex, respectively, and Ti has 
the largest dependence on temperature with an error value of 
o"A7"=0.17 dex. From the lower left panel of Fig. |4]we see that 
the errors associated to the uncertainties in metallicity are quite 
similar among the elements with cr^e/H] values ranging be- 
tween 0.05 and 0. 10 dex for an adopted 0. 1 dex error in [Fe/H]. 
Contrary to the sensitivity of the other elements, Mg is more af- 
fected by changes in metallicity than in temperature. The reason 
for such a dependence is most likely given by the blend with the 
strong Fel line at the red side of the triplet. Additionally Fig. [4] 
shows that the errors arising from uncertainties in gravity, o"i ogg , 
and microturbulence, cr^, are small. However, Si still shows a 
0.05 dex variation given an uncertainty of 0.3 dex in gravity. This 
is an important point to notice given that in our analysis we as- 
sume a distance of 8 kpc for all bulge stars, neglecting the bulge 
distance spread. If such stars are actually at 6 kpc, the bolomet- 
ric magnitude will be underestimated by 0.65 mag. Such an error 
will change the photometric gravity by 0.25 dex. Mg abundances 
are almost completely in-sensitive to changes in both log g and 
£ while Ca and Ti show a dependence of 0.06 dex on a microtur- 
bulence change of 0.2 dex. 

5.2. Errors from the spectral synthesis 

As described in Z08, stars in the field at b= -12 were observed 
twice due to a mistake in the fiber allocation. For this reason we 
only have about half the number of stars than in the other fields. 
However, we can use this repeated observation to carry out the 
same procedure to re-derive the abundances for all the elements 
analyzed in our work. In this way, we can estimate the influence 
of the steps in our analysis which are done by eye such as con- 
tinuum placement and FWHM adjustment used for the synthesis. 
Additionally this will also provide an estimate for the errors as- 
sociated to different S/N among stars in the same field. Figure [5] 
shows the difference between the two independent analysis. The 
mean value of the difference for all elements ranges between 
and 0.02 dex, which shows the absence of systematics in our 
procedure. The scatter is quite constant among the elements up 
to a value of 0.1 dex. Since the stellar parameters were not re- 
calculated, this scatter is not produced by their uncertainties but 
from the procedure to measure the abundances. We can therefore 
assume a 0. 1 dex error from this source, which added to the er- 
rors from the stellar parameters gives an estimation for the final 
errors in our measured abundances following the relation: 



[X/Fe] 



2 , _2 

+ °"A[Fe/H] + °"Alog£ 



(3) 



Therefore, the final values of the estimated error in our anal - 
ysis are 0. 15 dex for [Mg/Fe], 0.16 dex for [Ca/Fe], 0.22 dex for 
[Ti/Fe] and 0.14 dex for [Si/Fe]. 

5.3. The role of spectral resolution 

It is also necessary to consider the effect of using lower resolu- 
tion (R~22,500) spectra from GIRAFFE in comparison to the 
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Fig. 4. Abundance errors associated with the uncertainties in effective temperature (upper left), surface gravity (upper right), metal- 
licity (lower left) and microturbulence (lower right). A change of +200 K was applied to T e ff, +0.1 dex to [Fe/H], +0.3 dex to logg 
and +0.2 dex to Model atmosphere were then created and abundances were re-determined. Filled circles show the change in 
abundance as A[X/Fe] = [X/Fe]-[X/Fe] ± A when changes were added to the original parameters. Empty circles show the change in 
abundance when the respective value is subtracted. Dashed lines show the mean variations <x (+) and cr ( _) when the uncertainty in the 
stellar parameter is added and subtracted respectively. 



high resolution of the disk sample of R>45,000. In order to 
address this question we have re-calculated all the abundances 
using high resolution spectra for a subsample of stars in each 
field that was also observed with UVES(R~45,000). The abun- 
dances were calculated following the same procedure as for the 
GIRAFFE sample. 

Figure [6] shows the comparison between the abundances for 
each star in all analyzed elements. No systematic offsets are ob- 
served in our analysis from the change in resolution. This is par- 
ticularly important to notice since it could be a source of system- 
atic shifts in the trends when comparing bulge and disk. 

In order to further check the existence of any systematics we 
have also redetermined the abundances in the disk sample after 
degrading the resolution of the spectra to the one of GIRAFFE. 
Figure|7]confirms that no differences are observed from a change 
in resolution. Only Si shows a mean difference of 0.04 dex which 
is negligible within the errors on the abundance itself. 



Additionally we have re-measured the abundances in our cal- 
ibrator stars Arcturus, //-Leonis and the Sun, using the spec- 
tra with the original high resolution. The results were com- 
pared with our abundance measurements from convolved spectra 
(Tableland Fig. O and no differences were found. We therefore 
conclude that our comparison of bulge stars based on lower re so- 
lution to those of the disk is valid and no systematic shifts should 
be applied to our results. 

6. Results 

6.1. Thick and thin disk sample 

In order to characterize the bulge alpha abundances and compare 
them to other galactic components we have checked carefully the 
homogeneity of the procedure. Stellar parameters for the disk 
sample were calculated following the same procedure described 
in Z08. A comparison between the stellar parameters obtained 
using our procedure and the values presented in AB10 is shown 
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Fig. 5. Difference between the abundances measured using the 
same lines for stars observed twice in the field at b= -12. The 
dashed line shows the mean difference between both measure- 
ments and the dotted lines the lcr around the mean. 
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Fig. 6. A[a/Fe]=[a/Fe] UVES -[a[Fe] CIRAFFE for our bulge sample 
plotted as black filled circles. The dashed line indicates the mean 
value of the difference and the dotted lines show the lcr range 
around the mean difference for each element. 



in Fig. [8] We systematically find a larger effective temperature 
with a mean difference of 130 K between the two methods. Also, 
log g values are systematically larger in our analysis with respect 
to AB10. These differences between our studies show once more 
the importance of a homogeneous analysis in which parameters 
for all samples are on the same scale. 

Using our determination of the stellar parameters, we have 
measured the abundances for Mg, Ca, Ti and Si as described in 
section 3. Figure[9]shows the difference between the abundances 
for the disk sample derived in this work and the values presented 
in AB10. Mg abundances appear to be on the same scale how- 
ever abundances of Ca, Ti show larger differences. 



Fig. 7. A[a/Fe]=[a/Fe]HR-[a , /Fe]£R for the disk sample plotted as 
black filled circles where the low resolution (LR) spectra were 
obtained by degrading the high resolution (HR) to the resolution 
of our GIRAFFE sample. The dashed line indicates the mean 
value of the difference and the dotted lines show the lcr range 
for each element. 



AB10 uses the abundances of a subsample of thin disk stars, 
with metallicity near to the solar value, as internal zero-points. 
For the same stars we find larger Ca and Ti abundances as shown 
in Table [3] However we remain with our zero -points based on 
the solar abundances of lAsplund et al.l (120091) . Table [3] should 
be used to compare our results either with models scaled to so- 
lar abundances or other works on different scales. Within our 
log-gf and adopted zero points, the homogeneity of our analysis 
remains unchanged. 

6.2. Alpha elements across the bulge, thick and thin disk 

The final [cr/Fe] trends for Mg, Ca, Ti, and Si are shown in 
Figure [10] to [T3] and listed in Table |4] The general alpha en- 
hancement of the bulge stars is observed clearly in all four ele- 
ments. The alpha enhancement is also seen in our globular clus- 
ter members and their mean abundances do not differ from those 
of field stars at their given metallicity. The scatter in the observed 
trends is different: Mg abundances show much smaller scatter 
than Ca, Ti and Si abundances. In particular at low metallici- 
ties, where there is almost no dependence on [Fe/H], there is a 
~0.09 dex scatter observed for [Mg/Fe] and ~0.13 dex for the 
other elements. As shown in section 5, the dependence of the 
measured abundance on the stellar parameter uncertainties also 
varies among elements and the observed scatter in abundances 
is within the estimated errors. Therefore, the different amount of 
scatter is largely due to the uncertainties in the parameters. The 
similarity and overabundance in the bulge and thick disk stars 
appears to vary among the elements and therefore trends should 
be compared separately. 

Mg is the element which shows the best defined trend among 
the elements analyzed here. This element shows the least de- 
pendency on the stellar parameters and can be considered as 
the most reliable of the alpha elements for chemical enrich- 
ment studies as, together with oxygen, it is expected to be en- 
riched exclusively by SNII explosions. The mean abundance of 
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Table 4. Mg, Ca, Ti and Si abundances for stars in the four bulge 
fields. Metallicities are those presented in Z08 and listed here for 
reference. The full table is available in electronic form. 
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Fig. 8. Difference between the stellar parameters obtained in 
this work and the values presented in AB10 calculated as 
AX=X t hiswoik-XABio- The dashed line indicates the mean value 
of the difference and the dotted lines show the 1 cr range for each 
element. 
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Fig. 9. Difference between the [a/Fe] values for Mg, Ca, Ti and 
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cates the mean value of the difference and the dotted lines show 
the lcr range for each element. In particular, disk giants used 
by AB10 to calculate their internal zero points are shown as red 
filled stars. 
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Fig. 10. [Mg/Fe] abundances in 4 fields of the bulge shown as 
red filled circles. Bulge globular cluster members are shown as 
black filled squares. [Mg/Fe] abundances for the thick disk stars 
are shown as blue filled triangles and as empty black circles for 
the thin disk stars. 



the bulge at [Fe/H]< -0.5 is [Mg/Fe] = +0.31 which only dif- 
fers by 0.03 dex from the mean thick disk at the same metallic- 
ity range with the bulge being slightly more enhanced than the 
thick disk. Such a behavior remains in the intermediate metallic- 
ity range of our sample (-0.5 < [Fe/H] < -0.2) where the bulge 
shows a mean [Mg/Fe] of 0.27 and is still only 0.03 dex more 
enhanced that the thick disk. The differences we have found be- 
tween these populations are much smaller than those found by 
L07, in which the bulge is up to 0.2 dex more enhanced than 
the sample of thick disk dwarfs from the literature. Further com- 
parison at higher metallicities are difficult as separation between 
thin and thick disk becomes difficult. All elements show that, at 
the high metallicity range ([Fe/H]> 0) both the bulge and two of 
the high metallicity stars classified as thick disk by AB10 have 
lower ratios, as low as those of the thin disk. 

On the other hand, [Ca/Fe] trends show a larger scatter which 
complicates the comparison, particularly at high metallicities. 
However, at lower metallicities, where actually the thick disk is 
well defined, differences in [Ca/Fe] between the mean of both 
populations are higher showing the bulge 0.09 dex more en- 
hanced than the thick disk at [Fe/H]< -0.5 and 0.04 dex at 
higher metallicities. Worth notice as well is that [Ca/Fe] shows 
a shallower decline with metallicity than Mg in good agreement 
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Fig. 11. [Ca/Fe] abundances in 4 fields of the bulge. Symbols are 
as in Fig.fTOl 



Fig. 13. [Si/Fe] abundances in 4 fields of the bulge. Symbols are 
as in Fig.fTOl 
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Fig. 12. [Ti/Fe] abundances in 4 fields of the bulge. Symbols are 
as in Fig.fTUl 
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Fig. 14. [a/Fe] abundances in 4 fields of the bulge calculated as 
the average between Ca, Mg, Ti and Si abundances. Symbols are 
as in Fig.fTOl 



with model predictions where Ca contribution is also expected 
from SNIa explosions. 

Results for [Ti/Fe] are similar to those observed in Ca, show- 
ing trends with a shallower decline with metallicity but also a 
larger scatter. The mean abundance is [Ti/Fe] =0.40 and shows 
a difference with the thick disk of 0.10 dex from [Fe/H]<-0.5 
up to intermediate metallicities. Only above solar metallicities 
the mean Ti abundances of bulge and thick disk are the same. 
However these differences are within the scatter of the Ti abun- 
dances in the bulge which is consistent with the larger error 
dependence on temperature. Moreover, the [Ti/Fe] ratios in the 
thick disk are actually within the lower envelope of those of the 
bulge. 

Finally in the case of [Si/Fe], trends for the disk and the bulge 
are almost undistinguishable as in the case of Mg, with the bulge 
showing a mean [Si/Fe] =0.35 at [Fe/H]< -0.5 and differences 
lower than 0.03 dex with the thick disk at all metallicities. 



Figure [14] shows the average a element abundances. Trends 
for the bulge and both thin and thick disk are very well defined 
and we can also use them to compare the population properties 
between our different fields along the minor axis and with the 
field at a larger galactic longitude of 1=5°. A discussion regard- 
ing such a comparison is addressed in section 6.3. 

In order to generally compare abundances between bulge and 
disk, we have followed the procedure adopted in AB10 and fit 
two independent slopes for two different metallicity regimes. 
However, instead of fixing arbitrarily the point where the knee 
is produced to a [Fe/H] value of -0.5, as done in AB10, we have 
produced a linear fit to a metal-poor region between -1.2 and 
-0.5 dex and to a metal-rich region between -0.3 and 0.2 dex. 
These regions where selected in order to i) consider metallicity 
regions in which the trends are well defined and ii) avoid taking 
into account too metal-rich stars in which the uncertainties are 
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Fig. 15. The left panels show [a/Fe] trends as a function of [Fe/H] in 3 bulge fields located along the minor axis. Best fit trends are 
shown for both [Fe/H] ranges, a metal-poor between -1.2 and -0.5 dex and a metal-rich between -0.3 and 0.2 dex, as well as the 
location of the knee in all fields. The right panels show the [a/Fe] distribution for each field. 
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Fig. 16. Same as in Figure[T31for 2 bulge fields along the major axis. 



Table 5. Best fit parameters for linear (a+m[Fe/H]) trends for 
[<Mg,Ca,Ti,Si>/Fe] abundances in the bulge and thick disk. 
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larger. In this way, we are now able to provide a value for the 
metallicity at which the knee is actually produced, calculated as 
the intersection of the linear fits of both regions. The best fitting 
relations for both metallicity ranges in all the bulge fields and 
the thick disk are shown in Tableland plotted in Figures [TBI and 

As a general result, when considering the mean abundances 
of all the elements analyzed, differences between the bulge and 
the thick disk do not exceed 0.05 dex for all metallicities. Trends 
between both populations are therefore indistinguishable within 



the errors of our analysis. Additionally we see that the knee 
in the bulge in Baade's Window is located at a metallicity of 
[Fe/H]=-0.44 and differs by less than 0.1 dex among the other 
bulge fields. This knee is also present in the thick disk distri- 
bution at [Fe/H]=-0.51, only 0.07 dex different to the location 
for the observed bulge downtrend. Therefore, within the errors 
considered in our analysis as discussed in previous sections, no 
significant differences are observed between the bulge and the 
thick disk abundances for individual elements. The similarity is 
even stronger when all elements are compared in the mean. 



6.3. Is there an a -element gradient In the bulge? 

The existence of a metallicity gradient in the bulge is considered 
a strong indication towards the classical bulge formation sce- 
nario and is not expected to be present in boxy or pseudo-bulges. 
On the other hand, if Galactic bulge is a mixture of two (or more) 
components, it could be expected that these populations should 
have gone through different formation histories. A component 
formed via mergers would show higher alpha elements than a 
structure dynamically formed from disk material which would 
have enough time for SNIa to explode and decrease the [a/Fe] 
ratio up to Solar values. These components might have different 
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relative contributions across the bulge and therefore their content 
of alpha elements might help to distinguish between them. 

Figures[T5]and[T6]show that in the case of the Galactic bulge, 
the relations between [a/Fe] and metallicity are undistinguish- 
able among the different fields analyzed here. From the obtained 
relations for each field and in the two metallicity regimes that 
we have considered here, we can actually see no significant dif- 
ference between the fields in terms of the IMF and the timescale 
given by the knee at which SNIa contribution becomes impor- 
tant. Differences between the fields are not larger than 0.02 dex 
at all metallicity regimes and the knee is observed at the same 
metallicity within 0. 1 dex. This variation is in fact within our es- 
timation for the uncertainties in [Fe/H]. Such small differences 
are fairly consistent with the errors considered in our abundance 
analysis and a striking similarity can be concluded both along the 
minor axis (Fig. ITSb and the major axis (Fig. [Toll. Also shown in 
Fig- H3 and Fig. [T6j are the [a/Fe] distributions for each of the 
fields. The behavior of such distributions are consistent with the 
metallicity gradient in the bulge in which the population of stars 
with higher metallicities ([Fe/H]~0) and lower a element abun- 
dances ([a/Fe]~0) becomes smaller at higher distances from the 
galactic plane. 

6.4. Correlation with kinematics 

Coupling the information from radial velocities, proper motions 
and metallicity for the same bulge FLAMES sample that we 
present here. lBabusiaux et al.l (12010b demonstrate that metal-rich 
component shows disk/bar kinematics and a metal-poor com- 
ponent has kinematics of an old spheroid. The formation time- 
scale for the old (and metal-poor) spheroid component is much 
shorter, than that of the disk/bar component, and therefore it is 
expected that the former is alpha enhanced, while the latter has 
lower alpha element abundance. We investigate this by combin- 
ing our alpha element abundances with radial velocity informa- 
tion from Z08. Figure [T7]shows the radial velocity dispersion for 
different alpha element abundance bins. We observe that alpha- 
enhanced stars have a velocity dispersion of ~90 km/s, which 
is approximately constant along the minor axis, while the stars 
with solar or slightly sub-solar [a/Fe] abundances have larger 
velocity dispersion in the inner fields and decreasing velocity 
dispersion at higher latitude. In the b=-12° field the most alpha- 
poor component disappears altogether. The behaviour of alpha- 
enhanced stars mimics very closely the (lack of) trend shown by 
metal-poor component identified as having spheroid kinematics, 
while alpha-poor stars follow the trend of metal-rich componen t 
identified with disk/bar kinematics bv iBabusiaux et aD (|2010). 
While these trends are interesting, they need to be confirmed 
with a more complete coverage of the bulge. 

7. Discussion 

We have confirmed the chemical similarity between the alpha 
element abundances in bulge and thick disk within the errors as- 
sociated to our homogeneous abundance analysis. An overabun- 
dance in alpha elements of [a/Fe] =0.30 dex up to a metallicity 
of [Fe/H] — 0.5 is observed in both components as an indication 
of a short timescale of formation, with a chemical enrichment 
dominated by the contribution of core collapse supernovae. At 
this metallicity the observed downtrend in the bulge reflects the 
instance in which the Fe contribution from SNIa becomes impor- 
tant. In the bulge both the trend at low metallicities and the value 
at which the knee is produced do not vary along the minor axis 
for latitudes b<-4. This can be interpreted as a formation process 
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Fig. 17. Dispersion of radial velocity as a function of [a/Fe], in 
bins of 0.15 dex, for the three fields along the minor axis (filled 
circles and left axis scale). Overploted as a dashed line is the 
[a/Fe] distribution for each field, which scale is given on the 
right. Globular cluster members in the sample were removed for 
this particular analysis. 



of the dominating population of the bulge being fairly homoge- 
neous, sharing a similar IMF and formation timescale. Figure[T6l 
shows that these trends are indistinguishable also for two fields 
at different longitudes (1=1° and 1=5°) at a fixed latitude of b=- 
4°, an indication that the homogeneity is present along the major 
axis as well. 

The chemical similarity of bulge stars and those of the so- 
lar neighborhood thick disk, at least for metallicities in which 
thick and thin disk can be clearly separated, is also remarkably 
intriguing. Furthermore, iBensbv et al.l (1201 Ol) demonstrated for 
the first time that inner disk stars show both the dichotomy thin- 
thick disk and the chemical similarity with the bulge. Either 
the bulge and thick disk went through similar formation histo- 
ries or the process which formed the bulge is directly linked 
to that of the thick disk. Unfortunately the formation history 
of the thick disk is currently a matter of debate. Therefore, as 
long as a complete understanding of the formation mechanism of 
the thick disk is missing, a correct interpretation for the bulge- 
thick disk connection will be hard to achieve. Particularly inter- 
esting, in the context of our results, are the scenarios in which 
both bulge and thick disk sha re a fast and ear l y form ation such 
as the scenario dis cussed in Bourna ud et al.l (120091) (see also 
Imm eli et al.l 120041) . They show the properties obtained from 
simulations in which a thick disk and bulge are formed inter- 
nally by instabilities of gas-rich clumpy disks in young galaxies. 
Such clumpy and tur bulent rotating disks ar e in fact observed in 
galaxies around z~2 (Elmegr een et al.l2 008). This shared forma- 
tion mechanism could explain the observed chemical similarities 
between these components. 

The main result our observations provide is that all fields in 
the bulge follow the same [a/Fe] vs. [Fe/H] relation, and that 
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metal-poor stars correspond to a population enhanced in a el- 
ements which shows homogeneous kinematics. These observa- 
tional facts taken on their own are in accordance with the clas- 
sical bulge formation through the fast dissipative collapse. This 
is in contrast with the recent literature claim s of pure disk/bar 
origin of the Galactic bulge (She n et al.ll2010h . which would re- 
sult in shallower [a/Fe] decline, more similar to that of the thin 
disk. The most metal-rich part of the population, that shows a 
downturn in [a/Fe] has formed on a timescale long enough for 
SNIa to start contributing significantly with iron peak elements 
and shows kinematics corresponding to a disk/bar component 
according to Babusiaux et al. (2010). This component progres- 
sively dissapears with increasing distance from the plane. 

Altogether, or results provide observational evidence for the 
coexistance of a classical bulge, within a bar-like component. 
Dual component bulges have also been ob served in external 
galaxies in the work of iPeletier et al.l (120071) which shows an 
old elliptical-like component dominating the Bulge at higher dis- 
tances from the plane, while a disk-like component is observed 
in the cent ral regions. In particu lar, we highlight the recent re- 
sults from [Williams et all (1201 ll) . They show that the bulge of 
a pure disk galaxy, NGC 3390, rotates perfectly cilindrically, 
has chemical properties resembling those of the disk and shows 
no metallicity gradient. On the other hand, the boxy bulge of 
NGC 1381 shows more complicated rotation properties, is al- 
pha enhanced with respect to the disk and shows a metallicity 
gradient. These properties for NGC 1381 strikingly resemble 
those presented here, in a star by star basis, for the Milky Way 
bulge . In term of simulatio ns, the chemo-dynamical model of 
Saml and & Gerhardl (120031) also predicts a composite bulge: an 
early collapse component and a metal-rich, alpha-poor popula- 
tion associated to a bar. 

Clearly, to settle the issue of the origin of the bulge and to sit- 
uate it in the context of the formation of the other Galactic com- 
ponents, a detailed analysis of the kinematics and stellar abun- 
dances is necessary in a more complete coverage of the bulge 
and in the inner disk of the Galaxy. 

8. Conclusions 

We have analyzed the abundances of Mg,Ca,Si,Ti in four fields 
of the Galactic bulge. We have carried out a homogeneous anal- 
ysis comparing [a/Fe] ratios for a large sample of bulge stars 
(650) with those of giants belonging to the thin and thick disk in 
the Solar vicinity. We also compared bulge [a/Fe] distributions 
across the different fields in our sample. Our conclusions can be 
summarized as follows: 

1 . We have confirmed the chemical similarity regarding alpha 
element abundances between the bulge and the thick disk 
when a homogeneous analysis is done and using only giant 
stars spectra. Both populations are enhanced compared to the 
thin disk and show a downtrend in [a/Fe] starting at a metal- 
licity of [Fe/H] — 0.4. At higher metallicities, close to solar, 
bulge [a/Fe] ratios drop to values similar to those of the thin 
disk. 

2. The trends of [a/Fe] ratios are indistinguishable among our 
three fields along the bulge minor axis. Analysis on a field 
at 1=4°, b=5° hints for a homogeneity also present along the 
major axis. The knee showing the downtrend in [a/Fe] is ob- 
served in all fields at [Fe/H]~-0.4. A more complete mapping 
of different bulge regions is necessary to confirm the conclu- 
sions regarding the major axis lack of gradients as well to 
explore the [a/Fe] distribution for inner bulge fields (b>-4). 



3. The [a/Fe] distributions in our different fields show that, 
while the population of metal-rich stars and low [a/Fe] ratios 
observed at b=-4° seems to disappear at larger galactic lati- 
tudes, trends and distributions at low metallicities ([Fe/H]<- 
0.2), as well as velocity dispersion, remain unchanged at dif- 
ferent latitudes. Such results are consistent with recent evi- 
dences for a dual nature of the bulge. 

Our results demand future work to be based on a more com- 
plete coverage of the bulge. A better mapping of chemical and 
kinematical signatures, coupled with the complete morpholog- 
ical characterization to be obtained from a multi-epoch photo- 
metric survey such as the VISTA Variables in the Via Lactea 
(V V V: iMinniti et a"Dl2010l) . could provide the final step towards 
the observational characterization of the bulge. 
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